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Abstract. We study the dissipation mechanism of a stochastic particle sys- 
tem for the Burgers' equation. The velocity field of the viscous Burgers' and 
Navier-Stokes equations can be expressed as an expected value of a stochastic 
process based on noisy particle trajectories (Constantin, Iyer, Comm. Pure 
Appl. Math, 2008). In this paper we study a particle system for the viscous 
Burgers' equations using a Monte-Carlo version of the above; we consider N 
copies of the above stochastic flow, each driven by independent Wiener pro- 
cesses, and replace the expected value with i times the sum over these copies. 
A similar construction for the Navier-Stokes equations was studied by J. Mat- 
tingly and the first author (Nonlinearity, 2008). 

Surprisingly, for any finite N, the particle system for the Burgers' equations 
shocks almost surely in finite time. In contrast to the full expected value, the 
empirical mean JUi does not regularize the system enough to ensure a 
time global solution. To avoid these shocks, we consider a resetting procedure, 
which at first sight should have no regularizing effect at all. However, we 
prove that this procedure prevents the formation of shocks for any N ^ 2, 
and consequently as A' — ^ oo we get convergence to the solution of the viscous 
Burgers' equation on long time intervals. 



has been studied extensively from several different points of view. Here u > repre- 
sents the viscosity, making the equation dissipative in nature. The inviscid Burgers' 
equation, (equation with = 0) is studied as the basic example of a scalar 
conservation law (see e.g. [SlIZ])- The Burgers' equation is also linked to the KAM 
and Aubry-Mather theories [HlIIl]- It is the simplest PDE that models the Euler 
and the Navier-Stokes nonlinearity. As such, it has been extensively studied as the 
first step in understanding the two key unresolved issues in fluid mechanics: turbu- 
lence and regularity of the Navier-Stokes equations in three dimensions. In the first 
category the objective is to characterize the statistical properties of turbulence [6] . 
In the second category the objective is to understand the regularizing mechanism 
of dissipation [T1[T4]. This paper falls into the latter category: we study the reg- 
ularising mechanism of a particle system for the Burgers' equations, analogous to 
the particle system for the Navier-Stokes equations developed in [4l[TT|. 
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1. Introduction 



The viscous Burgers' equation 




dtu + udxU — vd'^u — 
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In ^ , a class of second order non- linear transport equations (the Navier-Stokes 
and viscous Burgers' in particular) were formulated as the average of a stochastic 
process along noisy particle trajectories. The formulation for the Navier-Stokes 
equations developed in j4j involves recovering the velocity u via the average of a 
non-local functional of the initial data. For the viscous Burgers' equation, however, 
the formulation is simpler. Explicitly, consider the stochastic flow 

(1.2) dXt = ut{Xt) + v^dWt 

with initial data Xo{a) = a for all a G M. Here W denotes a standard ID Wiener 
process. If we require that the velocity u satisfies 

(1.3) Ut = E [uo o (Xr')] 

where E denotes the expected value with respect to the Wiener measure, then u 
satisfief0 the viscous Burgers' equation (II. ip and initial data uq- We clarify that 
in (jl.3l) and subsequently, for any given time t ^ 0, denotes the spatial inverse 
of the diffeomorphism Xf. Namely, we know [18, Theorems 4.5.1, 4.6.5] that for 
regular drifts u, the stochastic flow X has a modification which is a stochastic flow 
of diffeomorphisms of M. Replacing X with this modification if necessary, for any 
t ^ 0, we define X^^ to be the inverse of the diffeomorphism Xf. That is, for any 
t ^ 0, we have Xt{X^^{x)) — x surely for all a; € E, and X^'^{Xt{a)) = a surely 
for aU a G R. 

Observe that when = 0, the system (|1.2p - (|1.3p is exactly the method of char- 
acteristics for the inviscid Burgers' equation. Indeed trajectories of the flow X 
are now characteristics, and equation (jl.3p states that the velocity is transported 
along characteristics. Thus, the v > Q case could be viewed as a stochastic gener- 
alization of the method of characteristics: we transport the initial data along noisy 
characteristics, and then average with respect to the Wiener measure. 

The usual Monte-Carlo method of solving (fr2|) - (IO)) numerically l20l[21] is to 
replace the flow Xt with N different copies Xl'^ , each driven by an independent 
Wiener process , and replace the expected value in p.3p by the empirical mean: 

J2f=i- Explicitly, the system in question becomes 

(1.4) dXi'^ = (^t'^) dt + V2^dWl 

(1.5) Xr(a)=a, 

(1.6) J^f = {Xl''')-\ 

(1-7) ^if^^oo^r 



where uq is the given initial data, a sequence of independent Wiener processes 
and v > the viscosity. As before, for any t ^ 0, {Xl'''^)~^ denotes the spatial 
inverse of X^' . Throughout this paper, with the exception of Section[3l we impose 
periodic boundary conditions on the above, and assume the initial data is periodic 
with period 1. 

For the Navier-Stokes equations, the particle system in |llj involves using a 
higher dimensional Wiener process, and replacing (jl.7p with the average of vorticity 



^This is only valid for spatially periodic or decay at infinity boundary conditions. 
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transport and the Biot-Savart law; 



(1.8) 



(1.9) 




i,N 



where = V x uq is the initial vorticity. In [11 , the authors considered the 
system (jl.4p - (jl.6p & (jl.8p - (|1.9p . with spatially periodic boundary conditions, and 
proved global existence in two dimensions, local existence in three dimensions, 
convergence to the correct limit as — > oo, and described the asymptotic behaviour 
for fixed iV as t — ?► oo. 

Surprisingly, the techniques of [TTj Jail for the particle system for the Burgers' 
equation (the system (|1.4p ~ (|1.7p ). Indeed, preliminary numerical simulations indi- 
cate that the system (|1.4p - (|1.7l) shocks almost surely, in time independent of N . 
We provide a class of initial data for which we can prove (|1.4p - (|1.7p shocks almost 
surely. We, however, we are unable to analytically prove that the shock time is 
independent of N . 

One heuristic explanation for the shock is as follows: This particle system (|1.4p - 
p.7p is dissipative only for short time jl_l, Theorem 5.2]. Once the system (|1.4I) - 
(jl.7p stops dissipating energy, the growth from the non-linear term should force the 
system to inherit properties of the inviscid Burgers' equation, which shocks if the 
initial data is not monotonically non-decreasing. 

We remark that the particle system for the Navier-Stokes equations (the sys- 
tem l|1.4p - (|1.6l) & (|1.8p - (ll.9p ) also dissipates energy only for short time 'TV, Figure 
1 and Theorem 5.2]. However, no dissipation is required to prove 2D global exis- 
tence for this system [TTl Theorem 3.5]. This is because (|1.8p - (|1.9p are structurally 
similar to Euler equations, for which 2D global existence is well known [2^ (see 
also i2jl9^). In contrast, however, the particle system ()1.4p - (|1.7p is structurally 
similar to the inviscid Burgers' equation which is known to shock in finite time. 

The natural approach one would expect to use 'overcoming' the shocks in (|1.4p - 
(|1.7p . would be to continue the system past shocks as weak solutions using an 
analogue of the Rankine-Hugoniot condition[7l Section 3.4.1], and then prove that 
as TV — )■ oo, these weak solutions converge to the smooth solutions of the Burgers' 
equation. This approach, however, is impossible to use as the stochastic PDE 
satisfied by involves second order terms, for which the classical techniques [71 
Section 3.4.1] will not work. 

While the system ()1.4p - (|1.7p can not be continued past shocks, the shocks can 
(surprisingly!) be "avoided with large probability" by resetting the Lagrangian 
maps. This is the main content of this paper. Namely, suppose we solve ()1.4p - ()1.7p 
for short time St, and then replace the initial data with u^, and restart the system 
(|1.4p - (|1.7p with this new initial data. Our main theorem shows that, if we repeat 
this procedure often enough, then we can avoid shocks on an arbitrarily large time 
interval, with probability arbitrarily close to 1. 

Explicitly, consider the system 



(1.10) 
(1.11) 




(1.12) 
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N 



(1.13) 



EN At. 



where A; £ N, and t is always assumed to be in the interval {kSt, {k + l)dt]- 

If St is small enough, we show that solutions to this system exist on arbitrarily 
large time intervals, with probability arbitrarily close to one. Once existence with 
large probability is established, it is easy to show that as — ?► oo these solutions 
converge to the smooth solutions of the viscous Burgers' equation. 

Before proceeding further, we remark that the fact that the shocks can be avoided 
by resetting is doubly unexpected! Firstly, we know that the inviscid Burgers' 
equations need to be regularized in order for them to have smooth solutions. The 
resetting procedure above should morally provide no regularization, as explained 
in Section[2] Secondly, the system (ll.2[) ~ (ll.3p is Markovian; if we reset it at regular 
intervals (as above), then the new solution obtained will be no different from the 
original solution without resetting. 

Fortunately, the system (|1.4p - (ll.7p is not Markovian, and if we reset often 
enough, the generic short time dissipative effect is strong enough to overcome the 
nonlinear growth, and with large probability prevents the formation of shocks. 
We observed numerically that even large resetting time 6t (i.e. comparable to half 
the shock time of the inviscid Burgers' system) is enough to ensure that the sys- 
tem (jl.l0|) ~ ()1.13|) is globally well posed. With the techniques in this paper however, 
we are only able to prove a global existence result for ()1.10|) - (jl.l3|) when St is small. 
The question for large St remains open, and cannot be addressed using techniques 
in this paper. 

Finally, we mention that our technique can be used to show global existence of 
the analogue of (|1.10p - (ll.l3p for the Navier-Stokes equations in two dimensions. 
As this is already known |11| . without resetting^ we do not carry out the details 
here. 

One interesting application would be to the 3-dimensional Navier-Stokes equa- 
tions. There are numerous results showing global existence of solutions to the 
Navier-Stokes equations with small initial data. One new, interesting question that 
can be asked in this framework is the existence of solutions for arbitrary initial 
data, which are time global for some small (non-zero) probability. The McKean- 
Vlasov type nonlinearity prevents us from asking this question for the stochastic 
Lagrangian formulation for the Navier-Stokes equations (4} equations (2.3)-(2.6)]. 
For the system (|1.4I) - (|1.6I) & (|1.8p - (|1.9l) . the empirical mean provides no 

regularisation, so it is unlikely to expect small probability time global solutions. 
However, the repeatedly reset version of (|1.4p - (|1.6p & (|1.8p ~ (ll.9p is free of the 
McKean-Vlasov nonlinearity, and is dissipative, making it a better candidate for 
small probability time global solutions. Unfortunately, there are obstructions in 
proving this result directly with the techniques used here, and we are working on 
addressing this issue. 

The plan of this paper is as follows: In Section [2] we establish our notational 
convention, and prove our main theorem. This proof relies on a few lemmas, the 
proofs of which we postpone to Sections HI [SI and [51 In Section [31 we provide 
an example showing that without resetting, the system (|1.4l) - (jl.7l) shocks almost 
surely. As mentioned earlier, once global existence is established the question of 
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convergence as — >■ oo is easily handled. We conclude the paper by studing this 
in Section [71 

2. The main theorem and its proof. 

Throughout this paper, we assume (il, E, P) is a probability space and use E to 
denote the expected value with respect to the probability measure P. Let iV ^ 2 be 
a natural number (which will be fixed throughout this paper), {J-t}t^o be a filtration 
satisfying the usual condition^ on fl, and W^, . . . , be N independent Wiener 
processes adapted to the filtration {JFt}t^o- We assume subsequently, without loss 
of generality, that = ^■ 

We use C'^(T), to denote the space of all periodic functions on M (with period 1) 
which have k continuous derivatives. We use Lp(T), H''{T) to be the Lebesgue p- 
space, and the Sobolev space of order s respectively, consisting of periodic functions. 
When writing norms of functions in these spaces, we drop T. For instance we use 
the notation ||u||_f/s to denote the H^{T) norm of u. 

We use a calligraphic script to denote the analogous spaces for processes, on 
random time intervals. Namely, given ^ 0, and a stopping time r such that 
T ^ Iq almost surely, we define 

C''([to,r];T) = {u \ u£ C° {[to,T];C^ {T)) a.s., and WrAt is J"t adapted} , 

C^{[tf),T]\T) ^ {u I u e C"([io, t]; i^(Tr)) a.s., and Ui-At is adapted} , 

-H^([to,T];T) ^ {u I u£ C°([to,T];iJ''(T)) a.s., and UrM is J't adapted} , 

where we use the abbreviation 'a.s.' for almost surely. For convenience, if r is any 
stopping time, not necessarily greater than or equal to we define C'^([io, t]; T) — 
C'=([to,r V to];T), and similarly for H^C^'. To avoid confusion with the C'=(T) 
norms, we explicitly use 

sup sup \\ut{uj)\\(jk 

to denote the C''{[to, r]; T) norm of u. 

We clarify, u E C''{[to,T];T) means that there exists an event ft' C ft with 
P(ri') = 1 such that Vw e il', i G [to, t{uj)], Ut{uj) G C'^(T) and wt(w) is continuous 
in t. Further, Vt ^ Ur/\t is J-( -measurable. In words, C'^{[to,T];T) is the set of 
all processes which have a C'^(T) valued, continuous paths modification and are 
defined on the random interval [to , t] . 

Note that our spaces involve processes which are continuous in time almost surely, 
and we are not interested in quantifying any further regularity with respect to time. 
When the regularity in time needs to be quantified, the definition the analogous 
spaces is not as elementary (see for instance ^15;). 

Our main theorem shows given any arbitrarily large T, we can make our resetting 
time 6t small enough so that a regular solution to (|1.10p - (ll.l3p exists up to time T 
with probability arbitrarily close to 1 . In order to formulate our Theorem precisely, 
we will need to define the notion of solutions to the reset system ()1.10p - (|1.13p with 
respect to a stopping time. This is our next definition. 

^By 'usual conditions' 1131 Definition 2.25] we mean that the filtration {J^t}t^o is right con- 
tinuous, and J-'o contains all P-nuU sets in J^oo ■ 
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Definition 2.1. Let to ^ 0, t be a spatially independent stopping time such that 
T ^ to almost surely, and utg be a C^(T) valued J^t^, -measurable random random 
variable. Suppose u € C^([io, ■'']; T) is a unique fixed point of the system 

/>rAt />rAt 

(2.1) XlfM^^+ Us{Xlf,)ds+ dWl 

Jto Jto 

(2-2) A^^,^(XI;^'' 
(2-3) = ]^II"to 

i=l 

Then we define 

For convenience, we adopt the convention that if r is any stopping time (not nec- 
essarily satisfying t ^ to), we define S^^'^utg = S^^'^^^^^utg. 

Remark. Note that it is essential to assume r does not depend on the spatial 
variable, as in this case if the drift u is spatially regular, then the process Xl'Jl 
admits a modification which is a stochastic flow of difFcomorphisms. Hence the 
spatial inverse ^^'^ = (X^'^)"^ is well defined. We will subsequently always assume 
our stopping times are spatially independent. 

Remark. In Lemma [2771 we will show that S^''!' is well defined. Namely, if fc ^ 1, 
and utg E C''{T) is J^tg-measurable, then Lemma 12.71 shows that there exists a 
(deterministic) time <i > to such that the process u defined by ut = Sf^^'l^Utg 
belongs to C'=([to, ^i]; T). 

Remark 2.2. Note that we can view the dependence of the operator 5'^^'.^ on the 
stopping time r as a dependence only through the time interval of definition. In- 
deed, for any fixed St and stopping time r, C^([0,t];T) solutions of (|2.ip ~ (|2.3p 
are unique up to indistinguishability. This follows immediately from a standard 
Gronwall type argument, and we therefore omit the proof. Strong uniqueness im- 
plies that the operator S^'^ satisfies a compatibility condition: for to ^ 0, con- 
sider two stopping times ti,T2 ^ to such that u} :~ S^^'J^Utg S C^([0, ti];T) and 
U2 ■= S^'l^utg G C^i[0, ^2]; T), then -u^ — before n A T2. That is has a modi- 
fication such that for all t ^ to, m^^^^^^^ = u]^^^/^^^. Thus, when the time interval 
of definition is clear, we sometimes omit the stopping time t as a superscript of our 
operator S. 

For notational convenience, we omit the first superscript N for the remainder 
of this section. Given a (spatially independent) stopping time t, a deterministic 
starting time to ^ 0, a C^(T) valued J^t^ -measurable initial data utg, and a resetting 
time St > small enough, we can define a C^([io, ''']; IT) solution of the (stopped) 
system (|1.10p - (|1.13l) iteratively by 

{Ut = Utg when t = to, 

4' = SI+^s.Am whenever te{to + kSt, to + (fc + l)St] for 
some k e NU {0}. 

We are now ready to state our main theorem. 
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Theorem 2.3. Let N > 1, T > 0, e > 0, s > 6 + ^, and suppos^ uq e i7-'(T). 
Then there exists St = ^t{T,s,s, \\uo\\h=), independent of N , such that for all St < 
5t, there exists a spatially independent stopping time r such that P{t > T) > 1 — e 
and the process u defined by (|2.4p (with to — 0) is in the space C®([0, t]; T). 

Remark 2.4. The compatibility condition in Remark 12.21 allows us to discuss the 
notion of a maximal stopping time Tmax for which the iterative procedure in (|2.4p 
is well defined. Consequently, Theorem 12.31 will show that this maximal stopping 
time Tmax is in fact at least T with probability at least 1 — e. 

We emphasize that the operator Stg^t is not a smoothing operator, which as men- 
tioned earlier is part of the reason why Theorem l2.3l is surprising. We can see St^^t 
is not smoothing from the fact that (|4.3I) , the stochastic partial differential equation 
(SPDE) satisfied by Ut = 5*0, tUo is not dissipative (15j . One can immediately verify 
this as the diffusive term in (|4.3|) does not necessarily dominate the noise. 

Another (perhaps more intuitive) way of understanding the regularity properties 
of Stg,t is via time splitting. The Stg can be time split into two parts: S}^, the non- 
linear solution operator associated with the inviscid Burgers' equations, and Sf^ 
the operator corresponding to resetting. By considering time split version of (|1.10l) , 
one can see that S^^ corresponds exactly to the operator 

i=i 

The operator 5^^^^ causes growth on the Fourier modes. It is well known that 
the damping provided by vd^, for any > 0, is enough to overcome this growth, 
and this gives us global existence for the viscous Burgers' equations for any strictly 
positive viscosity. Thus if the operator Sf^ provides damping comparable to vd^, 
then the usual methods can be used to prove Theorem 12.31 However, the operator 
Sf^ provides no damping, as can immediately be checked from (j2.5p : the operator 
norm of S^^ is exactly 1 (surely) in all Sobolev and Holder spaces. This is the main 
difficulty in proving Theorem 12.31 

We overcome this difficulty by considering the limit v := lims^^ou^* . It turns 
out that V satisfies a dissipative SPDE, and if the initial data is regular enough we 
obtain convergence in a strong norm of u''* to v. This is the key to the proof of 
Theorem 12.31 and is formulated below. 

Lemma 2.5 (Key Lemma). ie07V > 1, /3 e NU{0}, Tq > to > 0, and t be a (spa- 
tially independent) stopping time. Let Ut^ be a C'^^^iT) valued J- tg -measurable ran- 
dom variable, and it*' G C''+^([io, r]; T) be defined by 1^^. Let v e C''+'^([io, r]; T) 
be the solution of the SPDE 

(2.6) dvt + vtd^vt dt - -dlvt dt+^Y. ^^t' = 0' 



The theorem, and proof, remain unchanged if we instead assume that uo is a (T) valued, 
Jx)-™easurable bounded random variable. 

^The proof of this lemma never uses the assumption N > 1, and is valid even for N = 1. 
However, for A'^ = 1, Lemma 12.51 is vacuously true as assumptions 1)2. 7|l and 1)2. 8| l will never be 
satisfied for non-constant initial data. 
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with initial data v\t=to = utg, and spatially periodic boundary conditions. Let tq 
(t V to) A To, and U be a constant such thau 



(2.7) 
(2.8) 



sup 
sup 



.St 



^ U a.s. 



WvtWci+fi ^ U a.s. 



and let wf* = — UrAt- Then there exists a constant C = C{[3, U,Tq), indepen- 
dent of N , 5t and t, such that 

(2.9) 



sup E 



1/2 



Our main interest in this lemma will be for /3 = 2, as it will enable us to obtain 
a C^{[to,To];T) bound on u from a C^{[to,To];T) bound on v. A C^{T) bound is 
all that is needed to continue a solution locally, thus controling the (T) norm of 
u with large probability, independent of St, will prove our theorem. Since (j2.6|) is 
dissipative, uniform in time bounds of strong norms of v are readily obtained. 

Lemma 2.6. Let > 1, s G N, wo G iJ*(T). There exists a process v G 
'H*([0, oo); T) which is a solution to the SPDE ()2.6p with initial data uq and peri- 
odic boundary conditions. Further, there exists a constant Vg ~ Vs{s, ||uo||ff=) such 
that 



(2.10) 



sup 

t>0 



almost surely. 

We remark that ()2.10p is an almost sure bound on a strong norm of v. The 
reason we are able to obtain almost sure bounds is because if we "multpily by v 
and integrate by parts" (or more precisely, apply Ito's formula to HwfUla), we obtain 
an equation with no martingale part! This is carried out in detail in Section [5l 

Lemmas 12.51 and 12.61 will now allow uniform in time control of a strong norm of 
u^* . The only remaining ingredient is to obtain a C^(T) local existence result, and 
guarantee that the inequality (|2.7p is satisfied uniformly in St. 

Lemma 2.7. Suppose Ut„ is a C^{T) valued Ft„ -measurable random variable such 
that there exists a constant ^ such that \\ut„\\c'^ ^ almost surely. There 
exists To — To{Ui) > to and a process u^* G C^{[to, To]; T) such that u^* is a solution 
to ([231) with T = To. 

If further for some n CzN, utg is a C"(T) valued Ft„ -measurable random variable, 
and there exists a constant ^ such that \\utg\\c^ ^ almost surely, then u^* 
is C"([0, T]; T), and further there exists a constant Un = Un{U^, n), independent of 
N and St, such that 



(2.11) sup uf ^ Un a.s. 

for all St < To. 

^The assumptions H2.7| l- H2.8| l can be weakened slightly at the expense of a lengthier, more 
technical proof. The weakened assumptions however, still require more than /3 derivatives. While 
replacing I l2.7ll -l l2.8ll with a condition involving only /3 derivatives would be of sufficient interest to 
warrant a more technical proof, reducing 4 + /3to4 + /3 — e only obscures the heart of the matter. 
Since sufficient regularity on our initial data will guarantee I l2.7l l-l l2.8l l anyway, we assume they 
hold and avoid unnecessary technicalities. 
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Remark. The existence time Tq above only depends on a the (T) norm of the 
initial data. However, on the existence interval, any additional regularity of the 
initial data is preserved. 

We are now ready to prove the main theorem. (The Lemmas 12. 5[ 12. 6[ and 12.71 
will be proved in Sections SI [5] and [6] respectively.) 

Proof of Theorem \2.3[ Let (5t > be a small time, to be specified later, and let 
St £ (0,(5t) be arbitrary. Given a stopping time r, we define the operator ^^'^^j 

by 



St-r 



mSt^t 



■ ° ^{m+l)St,(m+2)St ° '^'m^t , (m+l)i5t ' 



where fc g N is such that kSt < t ^ (fc + l)6t. 

Let vt be the solution of ()2.6|) . By Lemma 12.61 and the Sobolev embedding 
theorem there is a constant Vi, such that 

sup \\vt\\ci ^ Vi 

almost surely. Let Tq = To(2Vi) be the local existence time in Lemma [^771 namely, 
for any initial data uq with ||wo||ci ^ 2Vi, and for any St < Tq, the process '.'^"mq 
is C^([0, Tq]; T). Without loss of generality we can assume that Tq is an integer 
multiple of St . 

Note that our assumption uq G Lemma [^751 and the Sobolev embedding 

theorem imply the assumption (|2.8p is valid for /3 — 2 (in this case, the supremum 
can in fact be taken over all t e R) . Similarly Lemma 12.71 guarantees that the 
assumption (|2.7p is valid for j3 ~ 2 and all St < Tq. Thus Lemma 12.51 can be 
applied. 

2Vi}. Then 



Let ill be the event {Hitji^Uci 



,p( 











< 









> 1 



^ 1 



St 



cs. 



1/2 



[Sobolev embedding] 
[Chebyshev's inequality] 

[Lemma 12. 5| , 



csT 

where the constant ci above is the constant arising in the Sobolev embedding 
theorem. An appropriate choice of St will make P{Vli) arbitrarily close to 1. We 
clarify that while our bound on P{^i) depends only on St, the event f2i depends 
on St- 

We define a stopping time ri by 



Ti(w) 



To if w ^ r^i, 

2To ifwef^i. 
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Note that by Remark E2 we have ^o^,V^'uo ^o*,V^"wo for all t e [0,To]. Thus, 
by the semi-group property, and the fact that To is an integer multiple of St, 

^s.,r, ^(S^of^'^o forte [0, To], 

° \^TlT°'^o:n'^o forte (To, 2To], 

as long as either side is defined. We claim that the right hand side above is well- 
defined and in C^([0, ti]; T). We see this as follows: First for t e [0,To], this is 
true by Lemma 12.71 and Remark 12.21 Now, for ut ^ ili and any t e [ro,2To], 
•^Tot^ '^^ j'^^f identity operator. Further for almost every a; e we have 
-^0!^°^^) = Uto^^) e C6(T) and ||m^^(w)||ci ^ 2Vi. Thus for almost any 
w e ni, and for every t e [To, 2To], J5^^;;[^m^^ (w) e C^iT) by LemmaO 

Using Sobolev embedding, Chebyshev's inequality and Lemma [231 as above, we 
can find an event 1^2 C fii such that P{^2) is arbitrarily close to P{ili). As before 
we define a stopping time T2 by 

^ti{lu) ifa;^f]2, 
3To if w e ri2- 



T2(w) 



and the solution u*' = S^qV^^uq e C^([0, T3]; T). A finite iteration will complete the 
proof. □ 

Finally we address the question of — > 00. For this purpose, we re-introduce 
the superscript of A^ to indicate the dependence on A^ of the process considered. 
Using techniques similar to [TT], we show that the solution of (12. 6p converges 
to the solution of the viscous Burgers' equation as A^ — ^ 00. 

Proposition 2.8. Let be the solution of p.6p with initial data uq, and u\ he 
the solution of the viscous Burgers ' equation ()l.ip with the same initial data. If 
Uq e H" , s > |, then for any T > 0, there exists a constant C = C{T,s, \\uo\\h^) 
such that 

„ II 1, Ari|2 C 

sup E\\u^-vt s; — . 

te[o,T] 

We prove proposition 12.81 in Section [T] We conclude by remarking that Propo- 
sition 12.81 Lemma 12.51 and an argument similar to the proof of Theorem 12.31 will 
show that for small enough St, as N ^ 00, u^'^* converges to the same limit on an 
event of almost full probability. 

3. Almost sure existence of shocks without resetting 

In this section we show that the system (jl.4p - (|1.7p develops shocks almost surely, 
for any A^. The existence of shocks is simpler to prove if we work with monotone 
functions on M, instead of periodic functions, and thus for this section only, we will 
work with (fr4)) - p7f|) on E instead of on T. 

Let T be a (spatially-independent) stopping time, and we interpret C^([0, t];R) 
solutions to ()1.4p - ()1.7p . in the natural way (analogous to (|2.ip - (|2.3p ). The main 
result of this section shows that even if we stop 'bad' realizations of (|1.4p - p.7p . 



we can never continue solutions past the time -rr-, — , unless we introduce a 

regularizing mechanism. 
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Proposition 3.1. Suppose uo G C^(R) is a decreasing function, and let u be a 
C^([0, t'];M) a solution of (|1.4|) - (|1.7|) with initial data uo- Then, almost surely, 

(3.1) ■ " 



T < 



||9xUo||j 



Remark 3.2. The numerically observed shock time, in the periodic case, is inde- 
pendent of N, and it is of the order 



||9xMol| 



with large probability. This indicates 

our bound (13. 1|) is far from optimal. 

Remark 3.3. One can sho\\|lthat as — >■ oo the solution to (|1.4l) - (|1.7p approaches 
the solution to (11.11) at a rate of However, it is well known that the solution 
to (|l.ip is smooth for all time and no shock develops, provided the initial data is 
for instance and bounded [7]. 

The numerics mentioned in Remark 13.21 however, indicate that no matter how 
large N is, the system (|1.4[) - (jl.7|) will only be a good approximation to the true 
solution of (11.11) for short time, in the order of tto — K . 

Remark 3.4. Monotonicity of the initial data uq is precisely the condition that con- 
strained us to work on the line instead of on the torus. Specifically, the assumption 
dxUo{x) < for arbitrary a: G K simplifies the proof of 13 . 1 1 considerably. Numerics, 
however, indicate that this monotonicity assumption is redundant, and (|1.4p - (|1.7p 
develops shocks for arbitrary (periodic) initial data. 



Proof of Proposition \3.1[ Assume for simplicity, and without loss of generality, that 
||9a;Uo||L°° = —dxUQ{0) = 1. Let the stopping time r be the first time t ^ t' such 
that d^Xl^'^iO) = 0. Explicitly, 



t' Ami{t I dxXl'^{0)^0} 



We will first show that, t ^ N, almost surely. 
Differentiating (|1.4|) in space gives 



(3.2) 



didxX^ 



1,N\ 



dxUf 



N 



dxXl'^ dt, 



for t < r', almost surely. Here, our notation dxU^ \-^i,n means 



dxU 



N 



X 



Xx)^dxunxl'''{x)). 



Differentiating equation (jl.7p in space, we obtain 



Q N 

OxUf 



X 



1 ^ 



Uo o A^' ) 



(3.3) 



N 

I—. 

i=l 



Uo 



Xt 



.OxA^ 



X 



for t < t' almost surely. Since by the chain rule. 



(3.4) 



dxA. 



1,N 



X 



^dxX^' = dx \A^' oXt' j = 1, 



^See for instance [111 Theorem 4.1], where the analogous result is proved for the Navier-Stolces 
equations. 
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multiplying equation p.Sp by dxX^ 

1 



l,Af 



(3.5) d^u 



N 



X 



N 



gives 

N 

dxUo + ^ dxUo 



4=2 



1,N 



for t < t' almost surely. 

Note that for a solution of the system (|1.4I) - (|1.7I) . for all i, the flow Xl 
M — )■ M is homotopic to the identity map via diffeomorphisms of R. The same is 
true for the inverse inverse A\'^ , and thus dx^t^ \x'^-'^(^xXl'^ > 0. Finally, since 
uq is assumed to be decreasing, we know that OxUq < 0, and thus equations (13. 2p 
and (P?5t yield 



for t < t' almost surely. This (ordinary) differential inequality, along with the fact 
that OxXq' = 1, necessitates t < N almost surely. 

Now, by definition of t, and continuity (in time) of dxX. 

(3.6) 



l,Af 



lim dxXl'^iO) = 



on the event {r < r'}. From (|3.5p and the chain rule we have 

dxUo 



dxUf 



N 



1 

N 



dxX, 



1,N 



N 

dxUo 

i=2 



oX. 



■dxJ^f 



for t < t' almost surely. Note that all the terms on the right have (|3.5p have the 
same sign. Thus if one of these terms approaches — oo, then necessarily the entire 
right hand side approaches — oo. Equation (j3.6l) immediately implies the first term 
approaches — cx) at a; = on the event {t < t'}. Hence, on this event we have 



lim 9,<(X^^(0)) 



oo, 



lim WdxUtWi^^ ^ 

t— >-! t- 

almost surely. Consequently, if w e C^([0, r']; R), we must have P(r < r') — 0. 
Hence t' = t < N almost surely. □ 



4. Proof of the key lemma (Lemma 12. 5p . 

In this section we prove convergence of u^' to v as St — s- 0. The basic idea is 
to show that the velocity in our reset system (|1.10|) - (jl.l3p satisfies the limiting 
SPDE (|2.6p with a small error which is controlled as <5t — >■ 0. 

By shifting time, we may assume without loss of generality that to ^ 0. Further 
replacing r with r A Tq if necessary, we may assume t = tq ^ Tq. Throughout this 
section, we adopt the convention that to = 0, and N, (3, Tq, t, tq, uq, m''*, v and U 
are as in the statement of Lemma 12.51 We also assume the processes X'^^.g^ . , A].^^ . 
are all as in (|1.10p - (|1.13p . and for notational convenience, we will omit the N and 
5t as superscripts throughout this section. 

We need a few lemmas before we can prove Lemma 12.51 In our first lemma we 
determine an SPDE satisfied by u on the interval (fc^t, [k + l)5t\. 

Lemma 4.1. We define the process to he the i*'' summand in (jl.l3p . Explicitly, 

i ]Uo fort = 0, 

Uf — < 

]wfc5t°^U,t forteiTAk6t,TAik + l)6t]. 



(4.1) 
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Then for a/H G {1, . . . , N}, the process G C''+''([0, t]; T), and satisfies the SPDE 

pTf\t 

(4-2) X{r^kS,} «At - "fe^t) + / {usd^cul ~ \dlul) ds+ 

JrAkSt 

+ / d^v^ dWl = 0. 

on the interval t G [kSt, {k + l)St]. Similarly, the process u £ C*^^([0, r]; T), and 
satisfies the SPDE 

(4.3) UrM - UrAkSt + / {usd^U,, - ds + 

■JrAkSt 




on the interval t £ [k6t, {k + l)<^t]- 

Remark 4.2. A more intuitive, though less precise, way of phrasing the SPDE's (|4.2p 
and (|4.3I) would be to say for t e {t Ak5t,T f\ {k + l)5t], u, satisfy the SPDE's 

dui + utd^u^ dt - \dlv\ dt + d^v\ dWl = 0, for alH G {1, . . . , N) 

1 ^ 

dut + utdxUt dt — \d^ut dt + — ^ dxu{ dWl — 
with initial data u\=kSt = UkSt and u\t=kSt = UkSf 

Proof. From |4',ITT] (see also [l6l[22]) we know that when r = oo, the process Aj.^^ . 
satisfies the SPDE 

dAls^^t + utdxA^kSt.t dt - \dlJ^k8t.t dt + dxAls,,t dWl = 

on the time interval (fc^t, (fc + l)5t]. Writing down an integral version of this in the 
presence of a stopping time, equations (|4.2p and (|4.3p follow immediately from (|2.3p 
and (|4.ip by a direct application of Ito's formula. 

To checlu u, w^, . . . , G C'^^^{\f), r]; T), note that continuity in time is immedi- 
ate. Further, the spatial regularity of u has already been assumed in the statement 
of Lemma 12.51 For u^, . . . , it^, note that the r and the noise are spatially indepen- 
dent in equation (|1.10p , and it immediately shows that each X^,^^ . (and hence each 
A\^ .) is as spatially regular as u, which in turn shows that each G C([0,t];T). 

□ 

Now we show that with a small error u satisfies the SPDE (|2.6I) stopped at r, 
and obtain bounds on this error. Let Ejr^^Y denote the conditional expectation of 
Y given FkSf Given any process /, and a stopping time t, we define the stopped 
increment A^/ by 

Afe/ = /TA(*:+l)5t ^ frAkSf 



The spatial regularity of u,u , ... ,u" follows directly from an assumption only on the initial 
data, and a standard iteration argument. This is contained in Section [B] However for Lemma l4.1l 
an iteration argument is unnecessary because of assumption l|2.7|l . 
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For the (deterministic) process /* = t, we define A^t by 

'5t ifT^ik+l)St, 
Alt ={TA{k + l)St) - (r A kSt) ^ Ir-kdt if kSt t < {k + l)dt, 

ifT<fc(5t. 
Finally, let L be the (non-linear) operator defined by 

Lu = udxU — ^d^u. 

Lemma 4.3. Suppose (|2.7I) holds for some /3 e NU {0}, and let e'/. be defined 



(4.4) 



/ 1 ^ \ 

4 = Alu + Luks, Alt + d^uks, E ^kW' j 

^ 7 = 1 ^ 



Then there exists a constant C = C(/3, ?7,To) (independent of N,k,dt and t) such 
that for all St ^ Tq and A: ^ ^ we have 

(4.5) supE\d!^e',ix)\' ^CSl 

(4.6) supE\E^,,^d^^s',ix)\^ ^CSl 

Remark. Since u and all derivatives of u are a priori uniformly bounded almost 
surely, the proof of this lemma is straightforward. Without this a priori bound, we 
would only obtain similar bounds on e'^^Jj^a and E\\Ejri^^^d!^e'f.gJ\^j^2, which 

are still sufficient for Lemma EH 



Proof of Lemma \4-3\ We assume throughout this section that C is a constant only 
depending on U and Tq which could change from line to line. Note first that 
assumption (j2.7p and equation (j2.ip immediately imply that for any id {!,..., N}^ 

sup \\dxXl\\ ^ C and sup 3_^^ < C 



almost surely. Now equation p.3p . immediately yields the same bound for u, inde- 
pendent of N. Thus, making U larger if necessary, we may assume without loss of 
generality that (|2.7I) holds for all the processes u, u^, u^, . . . , u". 

For any fc G N U {0}, t e {kSf, (fc + l)St], ni^ (3 + 2, differentiating (|L3)) n times 
gives 

and hence 

/ i'Tf\t \ 2 

\JTAkSt ' 



^In equation I I4.4I I. technically, Lu^-g^ is not defined when t < k5t. However, in this case, 
AJ.t = 0, so the value of Luf^s^ does not matter. We use this convention subsequently without 
further mention. 
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Note that 

JrAkSt J kSt 

Thus using (|2.7p for both u and u*, for any t G {kSt, (k + l)St] we have 

(4.7) supE\d'^UrM{x)-d^Ur^ksA^)\' i^C is^ + ^Y.^*] 

where as usual the constant C may change from hne to hne, provided it only depends 
on /3, U and Tq. 

Similarly, using (j4.2p and the above argument we have 

(4.8) supEx.ks,^^} |5"<At(a;) - d^UksAx)\^ ^ C^t- 

for any t € (k5t, {k + l)6t] and n < 2 + /3. 

Now, from the definition of e' and equation (j4.3p we have 

^rA(fe+l)5t -, W ^rA(fc+l)5f 

4 = - / - ^ E / ^-^i '^^i + 



\ J = l ^ 

^rA(fe+l)<5t -, N „TAik+l)5t 

(4.9) = / (Lufci, - Lm,) + _ V / - rfW^i 

JT/\kSt JT/\kSt 

almost surely. For the Ito integrals in the second term above, 

/ PTA{k + l)St \ 

E^.s, / {d.uks, - d^ui) dWl = 

\ JrAk&f J 

X{r^k5,}X{s^r} (dxUkS, " 9x<As) J = 0, 



/ TAkSt 

\Jfe<5t 

and hence 

rA(fc+l)>5t 



\ JrAkSt 

I E [X{kSt<ir}X{s<ir}^x {Luks, " LurAs)] ds 

J kSt 



{LukSt - Lus) ds 



kSt 

(k+l}St 2 

E [X{s^t}^x (LUrAkSt - LUrAs)] ds 

kSt 

where the last inequality follows from (14. 7p with n — 2 + /3. This proves 

For (|4.5p . note that the expected value of the square of the first term in (|4.9p 
has already been bounded by CS^ < CS^. For the second term, the Ito isometry 
gives 
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1 r{k+l)St 2 



and using (|4.8p with n = 1 + /? the proof is complete. □ 



We now prove that a time spht version of the SPDE \2.Q\ satisfies the same error 
estimates as in Lemma 14.31 



Lemma 4.4. Suppose (|2.8I) holds for some /3 G N U {0}, and let e'l he defined by 
(4.10) e'l = + Lvks, Al + d^v^s, ^ ^ A^M^^ 



T/ien t/ie bounds (|4.5p and (|4.6p /loW /or ej! . 
Proof. First note that 

?^rAt - Vr/\kSt = - / i^'s ~ 1^ / ^^^^ 

JrAkSt JrAkSt 

almost surely. Thus for any n ^ 2 + /? and t e [kSt, {k + l)(5t] using p.Sp gives 
(4.11) sup£; |5>,At(a;) - c)>,Afc5,(a;)|' ^ CJ*. 

Similar to the derivation of (|4.9p we obtain 

rA{k+l)5t -| /.rA(fe+l)5t 

(9x^^/c5t - l9a;Ws) dWl 



^fe = / (-^^Wfeit - Lvs) ds + —y] 

JrAkSt j^i JrAkSt 



from definition (j4.10p . The remainder of the proof is now identical to the proof of 
Lemma 14.31 □ 



We are now ready to prove Lemma [2.51 We remark that, the assumptions (12. 7p 
and (|2.8p are stronger than necessary. We only need 



(4.12) sup {\\ut\\ci+ii + Iktllci+z^) =^ a.s. 

(4.13) sup sup E (\dl+^ut{x)\^ + \dl+f'vtix)\') C/, 



and the bounds on e' , e" provided by Lemmas 14.31 and 14.41 above. The proof we 
provide below depends only on these weaker assumptions. 



Proof of Lemma \2.5\ Let ek = e'j, — e'^, where e'j,, e'J. are defined by Lemmas 
and 14.41 respectively. Using equations (|4.5I) , (|4.6p and the corresponding estimates 
for ej!, we have 

(4.14) suvEd^Ekixf i^CS^ 



-'t 

|2 



(4.15) sn^E\Er,,^d^,ek{x)\ ^ CS: 
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for all fc ^ yi-. As before, we assume C is a constant that only depends on /3, U and 
To, which may change from line to line. Now, estimates (|4.14|) and (|4.15p imply 



E\\ek\ 



(4.16) 

(4.17) ^ 11^^,,, 

For the remainder of the proof we will use the weaker estimates ()4.16p and (|4.17p . 
Now, recall Wt = WrAt = UrAt — VrAt, and we know wq = 0. Thus 



(4.18) = -af (LurAkS, - LvrAksJAlt ~ d^+'wks, E ^^^') + 

We first estimate E{d^Al.w)^ where k is any integer such that kSt ^ Tq. 

For this, independence of W^, the mean square of the matringale term in (|4.18p 
is bounded by 

2 



E 



/ 1 ^ \ / 1 ^ 



3 = 1 



N 



(4.19) 

Next, for the mean square of the first term in (j4.1 

> 2 



E (af (LUrAkStix) - LVrAkSt[x))Y 

^ CE [{dP,+^u,AkS,? + [dt^VrAkS,?] + 

+ CE [{dl.{UrAkStdxUrAkSt)Y + {.dl.{VrAkStdxVrAkSt)f 

^ Csup max E (\dl+^UrAksAx)\^ + l^l+'^WrAMt (a;) |^) -f 



C max 



(^WrAkStW^C^+fi + ll^'rAfe<5tllc'i+'3) 



Hence, using (|4.12p and (|4.13l) we obtain 

(4.20) E (af {LurAkS, - LvrAksJAltf ^ 6^C. 

By (|4.16p the mean square of the last term in (|4.18p is also bounded by CS^. Thus, 
squaring (|4.18p . taking expected values and using Young's inequality gives 



(4.21) 



3+1 



Wk5j^ + C6^t+3E (afefc) 



Now for any K s^^, 






K-l 










K-l 




= 2 ^ d^^wu, 


^ - {Lu 


fc=0 





K-l 



fc=0 



9f (LUrAkSt - LVrAkSt) - 



i=i 
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K-1 



+ ^(9fA^^.)^ 

A;=0 

Taking expected values, integrating in space using (|4.21|) and (|4.16p gives 

K-l „ 

(4.22) E\\dl^wKSi\\^2 ^ -2(5i V £; / d^WksMiur^kStd^rUrAkSt- 

-VrAkStdxVTAkSt) dx + 

K-1 

k=o 



fc=0 



wks,fdx + CK6f 



For the first term on the right of the inequality ()4.22p note 

(4.23) d^WkStd^ (UrAkStdxUrAkSt - VrAkStdxVrAkSt) = 

d^WkStd^ {WkStdxUrAkSt - VrAkStdxWkSt) ■ 

Observe that the mass (spatial mean) of solutions to equation (12. 6p is constant in 
time. The same is true for solutions to (|4.3p . Thus, for all t ^ Tq, J^u^Atdx = 
Jj uq dx — Jj VrAt dx, and hence /rj, Wt dx — 0. Thus integrating (|4.23l) in space and 
using the Poincare inequality, the term involving u above is bounded by 



dx'^kStd^ {wk&AxUrAkSt) dx 



^ C Wd'^WUtf^^ \\UrAk5t\\cf> + ^ 



For the term involving v in (|4.23|) . when all the derivatives fall on w we have 

d^WkSt^TAkStdx^^WkSt = ^VrAkStdx {df^Wksy , 
and if we integrate by parts, we can avoid the extra derivative on w. Thus 



d^WkStdx [vrAkStdxWkSt) dx 



^ C\\d!^Wk5t\\l2 hrAkStWcl^ 



Thus using (|4.12p and the above estimates, the first term on the right of (|4.22p is 

bounded by 

(4.24) 

- 26tE / d^WkStd^ {uTAkStdxUrAkSt - VrAkStdxVrAkSt) dx ^ CdtE \\d^WkSt\\1^2 ■ 

For the second term in (|4.22p . we know WkSt is J^fe^j -measurable. Thus us- 
ing (j4.17p and the Cauchy-Schwartz inequality we obtain 



(4.25) 



/ \Edl^WksMek\ dx = f \E{d^^Wks,E^,s,d^,ek)\dx^CS'/\ 



The third term on the right of (|4.22p is always non-positive, and can be ignored. 
Thus, recalling K ^ and using equations (|4.24p and (|4.25p in equation (I4.22p 
we have 

K-1 

WkSt T2 ot- 

fc=0 
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The remainder of the proof is an elementary discrete Gronwall argument. Let 

K-l 

Then 

Vk+i - Vk = CStE \\d^WkSt 11^2 < CStUk 

and hence 



1/2 

Iterating, and using yo — C5^ gives 



Since fc ^ ^ this gives 



j/fc+i =^ (1 + C5t)vk- 
gives 

Vk ^ {l + C6tfC5l'\ 



(4.26) max yk ^ CS]'^ sup (1 + CS'^f''/^'* (7(5^^6^^° . 

5;>o 

This proves ()2.9p for all times t which are an integer multiple of 5t. Since for any 
a; S T, and fc ^ ^ we elementarily have 

sup E \ d^VrM{x) - ^'^Vrr^k5t{x)\^ ^ C5t 

kSt^t^(k+l)5t 

and sup E \d'^UrM{x) - 9f WrAfcSt (a;) f ^ C5t 

kSt^t^(k+l)5t 

completing the proof. □ 

5. Proof of Lemma [2T6] 



In this section we establish uniform in time bounds for the solution of (|2.6p and 
prove as in Lemma 12.61 We do this via the following two Lemmas: 

Lemma 5.1. Let uq g C°°{T), T > 0, and suppose v G C°°{[0,T];T) is a solution 
to (j2.6p with initial data uq and periodic boundary conditions. Then for any s G Z"*", 
there exists a constant Vs = Vs{s,T, \\uo\\h^) such that 



sup WviWhs S% Vs 



almost surely. 



Lemma 5.2. Let uq e C°°(T), and suppose v £ C°°([0, oo), T) is a solution to (|2.6p 
with initial data uq and periodic boundary conditions. Then for any s G T > 0, 
there exists a constant Vg — Vs{s,T, \\uo\\l2) such that 

(5.1) sup||i;t||^. sc:^ 

almost surely. 

We draw attention to the fact that apriori bounds are almost sure! Indeed, 
applying Ito's formula to ||wt||/,2 immediately yields an equation with no martingale 
part (see equation (|5.2p below). 

Given Lemmas 15.11 and 15. 2[ the proof of 12.61 is now immediate. 
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Proof of Lemma \2.6i Given the almost sure apriori bounds in Lemnias l5.1l and l5.21 
existence of solutions to ()2.6p follows via standard methods. The time global 
bound (|2.10p is also an immediate consequence of Lemmas 15.11 and 15.21 □ 



We devote the remainder of this section to proving Lemmas 15.11 and 15.21 



Proof of Lemma \5.1\ We prove Lemma 15.11 via energy estimates. First note that 
Ito's formula and (|2.6|) give 

1 ^ 

d{vtf = 2vtdvt + j;r2^{dxvtf dt 



-2vld^v^ dt + vtdlvt dt - 2^^ + J^^^-^tY dt. 



Integrating in space, and using VtdxVt dx — — v^dxVt dx gives 



(5.2) dt\\vt\\l. = -^--j\\dxVtrL^ =^ WvtWL^^Wu^L^ 
almost surely. 

A similar calculation shows ||wt||LP ^ I|wo||lp for all p ^ 2, and henc^ |lwt||L°° ^ 
IImqIIloo. Recall s > 1 by assumption, and so the Sobolev embedding theorem shows 
||wo||l°° =^ c||Mo||_ff for some absolute constant c. 

Now, differentiating (|2.6I) with respect to x and applying Ito's formula to {d^VtY 
we obtain 



1 



d{dxVtY = 2dxVt d{dxvt) + - \dlvt\ dt 

-2dxvt I dx (vtdxVt) dt - dt+^ 



j2dW^^+j^\dlv,\' dt 



Integrating with respect to x on [0, 1], and noting that dxVtdxVt dx — 0, gives 



d\\dxVt\\\2 = ^ ( ^ ^ ) ll^^'^tllf.^ dt+[2 j d'xVtivtdxVt) dx ] dt 



1 



1 



dt \\dxVt\\l2 ^ ||9^f;t||^2 + 8 \\vtdxVt\\l2 



(5.3) ll^^wtll^, + 8 lluollioc, \\dxvt\\l2 , 
almost surely. Thus, (|5.2p . (|5.3p and Gronwall's inequality gives 

(5.4) \\vt\\Hi^Cie'°' and f \\vt'\\l2 dt' ^ Cie'°\ 

Jo 

almost surely, for some constants Ci = Ci(||uo||ffi) and cq = co(||uo||l°°). 

For the remainder of this proof we adopt the convention that c, C denote absolute 
constants, Cs — Cs{s, ||uo||h=) denotes a constant depending only on s, ||uo||ffs and 



^This can alternately be shown using a version of the maximum principle | 17| . 
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Co denotes a constant depending only on ||uo||loo . The exact value of these constants 
are immaterial, and we will allow them to change from line to line. 

Similar to (|5.3|) . differentiating (12.61) twice with respect to x, applying Ito's for- 
mula to {dlvtY, integrating in space, noting /^9^ut9|wt dx = Q and using Holder's 
inequality gives 



dt Ptvth, ^ _ 1 _ ptvth, + 2 \\dtvth, \\d. {vtd.vt)\\^. 



^ -c\\dlvt\\\, + c{\\ dxVt\\\^ + IK'tllioc ) \\d; 



2 ||2 



^~c\\dlvt\\U + c\\vt\\lAdlvt\ 



almost surely, where the last inequality is obtained by the Sobolev embedding 
theorem. Using ()5.4p . this gives 

WvtWn^ < C2e^oll''*f„2<it ^ C-seC-i^^"' and /* \\vt,\\l, dt' ^ C^e''-^'°\ 

JO 

almost surely. Proceeding inductively, suppose we know 

\\vt\\H^ ^ exp (Cs-i exp {Cs-2 ■ ■ ■ exp {cot) ■ ■ ■)) 



(5.5) 



ft 

ut'||^,+i dt' < Cs exp [Cs-i exp (Cs-2 ■ ■ • exp {cot) ■■■)) 



holds almost surely for some s e Z+ . Differentiating (|2.6I) s + 1 times with respect 
to a:, applying Ito's formula for (9|+^wt)^, and integrating in space we obtain 

= - (i-^) Px-^'^t^L^ dt+2\\d^^\t\\^Ad:{vtd,v,)\\^2 dt 

since ^jdp'^vtdp''^Vtdx — Q. Thus 

dt pl+^vtf^, ^ -c \\dl+^vt\\\2 + C(\\dlvt\\l^ WdxVtWl^ + ■■■ + 

+ \\dxvt\\l^\\divt\\l2 + \\vt\\l^ ll^^'^*!!!^) 

< -c \\di+^vt\\\2 + c(||a>i||i^ + • • • + Wd^vtWl^ + 

almost surely. Thus by Gronwall's Lemma 

^ Cs+iexp(^J^ Wvt'fns+i dt'^ 

^ Cs+i exp (Cs exp (Cs-i • • ■ exp (cot) ■ ■ ■ )) 

almost surely. Further 

/ ||wt'llH=+2 rfi' < Cs+i cxp(Csexp(Cs-i • • -exp (coi) • • • )) 7 
Jo 

almost surely, completing the inductive step. By induction, (|5.5p holds for all 
s e Z+ completing the proof. □ 
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Proof of Lemma \5.2l We prove Lemma [5.2l via a bootstrapping argument in Fourier 
space. To fix notation, for n E Z, we use f{n) — e^^'^™^/(x) dx to denote the 
n*^ Fourier coefficient of /. 

On Fourier coefficients, using ud^u = ^d^u'^, equation (j2.6l) reduces to 

2nin ^ 

(5.6) dvt{n) H — j^Vt{n) dW^ + 2Tr^n^Vt{n) dt + win Vt{n — m)vt{m) dt = 

j=l mGZ 

for every n eTj. 

By Ito's formula applied to (|5.6p 



d|-Ot(n)| = vt{n) dvt{n) dt + vt{n) dvtin) -i |wt(n)| 

(5.7) = -4^2„2 (^1- |f;t(n)|' + mn (vt{n)Btin) ~ iAnjBtin)^ dt, 

where ?)t(n) denotes the complex conjugate of Vt{n), and 

Bt{n) = ^ {}t(n - m)vt{Tn), 

is the non-linear Fourier coupling in (15. 6p . Using > 1 and Young's inequality 
in (|5.7p gives 

dt \vt{n)\^ ^ -2TT^n^ \vt{n)\^ + 2TTn \vt{n)\ \Bt{n)\ 

(5.8) ^ -cn^ \vt{n)\'' + C \Bt{n)\^ 

almost surely, where, as before c, C are absolute constants (independent of uo,T), 
which may change from line to line. Thus, for any tg^O we have 

ft 

(5.9) |{)t(n)|' s; |t;t-(n)|'e-"'=* + C / e"™''*"*') |Bt,(n)|' dt' 

Jt'o 

almost surely, by Gronwall's inequality. 

By Parseval's identity we know (n)| ^ ||'yt|li2, and by conservation of energy 
(equation (|5.2p ') this gives \Bt{n)\ < ||uo|li2 almost surely. Thus the second term 
in the previous inequality is bounded from above by ^^^111*011^2. Since |?it;,P ^ 
||tt(/^||^2 ^ ||uo||^2, given a lower bound on t — t^, we can certainly arrange the same 
inequality for the first term. Thus choosing ti = for instance, and applying (15. 9p 
with = 0, we obtain 

(5.10) suv\vt{n)\' 

almost surely, where Co — Co(|lwo||i2, T) is some constant. 

Now we bootstrap, and use (|5.10l) to obtain a better estimate on Bt. Assume 
inductively that for some a G Z+, and ta = ■^;^T, we have 

(5-11) sup \vtin)f < -rkrr 

t^ta \n\ 

almost surely. Here Ca — C'Q(|juo||L2,T, a) is a constant which we allow to change 
from line to line if necessary. We will now establish (|5.1ip for a + 1. Note that 
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almost surely, for any t > ta, we have 

\Bt{n)\ \vt{n - m)\\vt{m)\ 1^ 2 ^ \vt{n - m)\ \vt{m)\ 



(5.12) 




J2 ^2\\uo\\ 

\m\^\n\/2 





1/2 



Now returning to (15.91) and choosing t'^ — t^, we see that the second term is 
bounded by i — ^a+i ■ For any t ^ ia+i, we can certainly arrange the same 

inequality for the first term, and hence this establishes (|5.1ip for a + 1. 

Finally note that if (|5.1ip holds for a, then ()5.ip holds for any s < f — 1, finishing 
the proof. □ 

6. Proof of Lemma [^771 

In this section, we prove the almost sure C" (T) bounds on u stated in Lemma [277l 
We need a few preliminary results first. 

Proposition 6.1 (Local existence without resetting). Let utg be a C^(T) valued 
J- to -measurable random variable such that 

almost surely. There exists Tq = Tq{Ui), independent of N , such that the solution 
to (|2.1|) - (|2.3|) exists on the interval [^Oj^o +?o]- Further if for some n ^ 1, Ut„ is 
a C"(T) valued, J- tg -measurable random variable with 

\\ut„\\c^^U° 

almost surely, then there exists U„ — Un{U^,n) such that 

(6.1) sup llwtllc" < Un 

almost surely. 

Proposition l6 . 1 1 can be proved using a standard Picard's iteration. A proof of the 
analogous result for the Navier-Stokes equations appeared in the appendix of [11] 
(see also [911101). The proof of 16. II is very similar, and we do not provide it here. 

Lemma 6.2. Let / ; R — > M denote the identity function, d G [0,1), and let 
A S C"(T) be a periodic function such that \\dx\\\L°° ^ d. Then there exists a 
constant c„_i — Cn-i{\\d^^^\\\L°° ,d,n) such that for any f G C"(R), 

(6.2) Wd^ [/ o (/ + A)]||^^ ^ liar/IL., (1 + ||5,A||^^)" + c„_i \\d:\\\^^ , 

(6.3) ||a^(/ + A)-i||^^ ^c„_i||a:A||^3., 

for n > 1. 

Remark. Note that since ||9a;A||Loo < 1, the function / + A is a C^(M) diffeomor- 
phism of R. The notation (J + A)~^ in (j6.3p refers to the inverse of the C^(R) 
diffeomorphisni / + A. 
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Proof. First note that we can view A as a periodic function (with period 1) in 
C"(K). Further, by the mean value theorem, for any fc ^ 1 there exists x G T such 
that a^A(x) = 0. Thus for any fc e {1, . . . , n}, we have |5^A| < c(n)||9;'A||Loo, for 
some constant c(n) depending only on n. (For fc = 0, we need to subtract the mean 
of A for this bound to be valid.) 

Now for any two f,g£ C"(M), we have 

n m 

(6.4) d:{fog) = J2{dTf)og ^ Hd^g. 

m=l ki-\ |-fcm=n»=l 

To prove (|6.2p . we set g = I + \. The term in (|6.4p corresponding to m — n gives 
the first term of 16.21 When m < n, we notice that ki > 1 for at least one i, and 
kj i^n~l for aU other j. Thus (/ + A)||l~ = ||9^'A||loo <; c(n)||a;^A||i=o . The 

remaining terms dx^ (I + X), j ^ i in the product can be bounded by c„_i. This 
proves (16.21) . 

For (ESI), set X = / + A and A = X~^. Since n > 1, a^(A o X) = 0, and 
using (|6.4p we obtain 

n — 1 m 



X {dxxy 



X 

feiH hfe™ 

/Ci>l 



By induction, one can assume that ||(9™A||loo ^ c„_i for all m ^ n— 1. Since d < 1, 
^ ^'^'^ remaining terms can be bounded by the same argument as 
before. This proves (|6.3p . □ 

Lemma 6.3. Let n € N, utg be a bounded, C"(T) valued, J- to -measurable random 
variable. For k € {0, . . . , k}, let be a constant such that \\uta\\ci' ^ almost 
surely. Let u be the solution of (|2.ip - (|2.3p with initial data ut = utg when t = to- 
If n > 1, there exists fl' £ with P{i^') ~ 1, To = To(C^i) > o.^'d- o, constant 
Cn-i — Cn-i{U^_i,n) such that 

(6.5) ||a^'ut(c^')llL^ U° (1 + c„_i(i - to)) 

for all Lo' e fl' , t e [toi^o + To]. For n = 1, (|6.5p holds with co to be an absolute 
constant. 

Proof. For simplicity, we assume to — 0. One can check that this assumption does 
not affect our proof below. Our first step is to obtain almost sure (T) estimates 
on the Eulerian and Lagrangian displacements. Throughout this section, we use 
the convention that c„_i = c„_i(C/°_]^, n) is a constant depending only on n and 
U^-i (or an absolute constant for n = 1), which can change from line to line. 

Let To = To{Ui) be the local existence time given by Proposition 16.11 and 
ci — ci{Ui) the almost sure bound on ||ut||ci from (|6.ip . Let / : M — > R be 
the identity map, X*, respectively be as in (|2.ip . ()2.2p . with t = Tq. Define 
XI = XI -I,£l^Ai- I. 

Differentiating (|2.ip with respect to x we obtain 

t 



Q.KWl^ ^ / Wd.UsW^^ (1+ ||9.A 
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almost surely, for t e [0,To]. By Gronwall's lemma, 

||9a;At||^^ e'^i* / ll^^u^ll^^ ds a.s. 
Jo 

for t e [0, To]. Recall t ^ Tq, ci only depends on Uq, and for all s ^ Tq, ||9a;Ms|lL°° ^ 
almost surely. Thus, as Tq is allowed to depend on [/°, by making Tq smaller if 
necessary we can arrange 

II 11 /■* II 11 1 

(6.6) \\dxK\\ro.^co \\dxUs\\^^ ds and sup p^^AJ < - 

Jo O^t^To ^ 

almost surely, for some absolute constant cq. Now 

1 (9:rAj) o 

" * ^ ^ " (a.X*) oAl^^^'l + (dxXl) o 
almost surely. Thus we must have 

(6.7) \\dJl\\^^^2\\dxXl\\^^ 
almost surely for t e [0,To]. Using (12. 3p and (|6.7p we have 

1 ^ 

Wd^^^th^ ^ ^ E ll^-"olLoo (1 + \\dJl\\^^) 
1=1 

1 ^ 

^]^EII^-"oIIl- (i + 2||a.Aj||^^) 

'i=l 

5% II<9:eUo||^^ +2co / WdxU^W^^ ds 
Jo 

almost surely for t G [0,To]. This proves (|6.5p for n = 1. 

For n > 1, local existence CProposition 16.11) guarantees that ||ut||c"-i ^ c„_i 
almost surely for t G [OjToli where c„_i = c„_i (J7^_j^, n). Assume by induction 
that the bound (|6.5p holds for some integer n — 1. This bound and equation (|2.ip 
immediately imply that ||9a;AJ||c»-2 ^ c„_i almost surelj|3 for t e [0,To]. Equa- 
tions (|6.6p and (|6.3p will imply ||9a;^t||c7i-2 ^ c„_i almost surely for t G [OjTo]- 

Thus using equations (|2.ip and (|6.2p we obtain 

lla^Ajll^^ ^ /V"K°(^ + ^s)]|Il~ 
Jo 

s=:c„_i r[||9>.|L=. + ||9;%|L^] ds 



almost surely. Using Gronwall's lemma this implies 

(6.8) II^^Ajll^^ s^c„_i fwd^u^h^ ds 



"'^"We remark that our somewhat unusual notation ||9a;AJ||pTi-2 instead of ||AJ||(-.n-i is neces- 
sary. This is because it is impossible to obtain almost sure bounds on ||AJ||jjOo . However, as our 
argument shows, we can obtain almost sure bounds on ||9^AJ||ioo for any fc ^ 1. 
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almost surely. Here we absorbed the constant e'^"-^* into c„_i, which is valid as 
< ^ To = Tq{U^). Now 



\\d>t\ 



1 ^ 
N ^ 

i=l 
N 



{\\d>o\ 



Id r 



^i\\d:ei\ 



1 



i=l 



< ||a>o|lico(l + c„_it)+c„- 



ds 



almost surely, where we used (|6.3p and (j6.7l) to obtain the second inequality, and 
equations (|6.6p and (|6.8I) to obtain the third inequality. Now Gronwall's Lemma 
gives (|6.5p . where we again absorb the exponential factor e'^"-^* into (1 + c„_ii), 
by replacing Cn-i with a larger constant, which by our convention we still denote 
by Cn-i- □ 

Proof of Lemma \S71\ By Proposition l6.ll existence will follow if we establish (j2.1ip 
for n = 1. We prove (j2.11[) by induction. Since the constant cq in Lemma 16.31 is 
absolute, the proof for n = 1 is identical to the proof of the inductive step. Thus 
we only prove the inductive step. 

Assume that (12. lip holds for n— 1, choose c„_i = c„_i([/„_i) to be the constant 
from Lemma 16.31 Thus whenever 5t <Tq, 



(6.9) 



"(fc+l)5t 



^ (l + c„_i(5t) 



a.s. 



holds for all fc ^ ■j'-. Iterating this we have 



for all t ^ Tq. Thus we choose C/„ to be given by 

Un = ||9"uo|Il~ sup (1 + C„_i(5) 
(5>0 



To/S 



From (|4.3p we see that u*' is conserved almost surely. Since uf* is periodic. 



bound on will give us a bound on ||it(' Hen, concluding the proof. □ 

7. Proof of Proposition! 



In this section we prove Proposition l2.8l We reintroduce an iV as a superscript to 
explicitly keep track of the dependence of our processes on iV, and prove convergence 
as iV — )• oo. 



Proof of Proposition \2.8[ Let w 



Then (HI]) and dH]) give 



(7.1) dwf + w^d^v^ dt + u\d^w^ dt - ^S^wf dt + 



N 

x Uj \ " 



dxv. 



N 



N ^ 



dWi = 0. 



Thus, by Ito's formula 



1 



d < %+ / [w'^Yd^v'^dx] dt+ 



■^w^ dxW^ dx ) dt-\-— pxW. 



N\ 
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Taking expectations and integrating by parts we obtain 



dtE\\w^\\"+E 



{w^f {2d^v^ - d^u\) dx 



By Lemma 12.61 and the Sobolev embedding theorem, there exists a constant C — 
C{s, \\uq\\h=), independent of N, such that 

sup ||9:cWt|lioo ^ C 

almost surely. It is well known that the same estimate holds for dxU^- Further, 
since £J||9j;Wf ||^2 supQ||9a;w|||oo , making C larger if necessary we have 

snp E \\dxVt\\l2 < C. 



Thus 



d,E\\w^\\l^v(E\\w^\\l, + l^ 



and, since wq = 0, Gronwall's lemma gives 

finishing the proof. □ 
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